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The radio light curves of SN 1993J are discussed. We find that a fit to the 
individual spectra by a synchrotron spectrum, suppressed by external free-free 
absorption and synchrotron self-absorption gives a superior fit to models based 
on pure free-free absorption. A standard r~ 2 circumstellar medium is assumed, 



and found to be adequate. From the flux and cut-off wavelength, the magnetic 



CO 

o 

field in the synchrotron emitting region behind the shock is determined to 
B pa 64 (i? s /10 15 cm) -1 G. The strength of the field argues strongly for turbulent 
amplification behind the shock. The ratio of the magnetic and thermal energy 
density behind the shock is ~ 0.14. Synchrotron losses dominate the cooling of 
the electrons, while inverse Compton losses due to photospheric photons are less 



O 

important. For most of the time also Coulomb cooling affects the spectrum. A 



Oh 
6 

model, where a constant fraction of the shocked, thermal electrons are injected 
and accelerated, and subsequently lose their energy due to synchrotron losses, 
reproduces the observed evolution of the flux and number of relativistic electrons 



well. The injected electron spectrum has dn/d'j oc 7 -2 ' 1 , consistent with diffusive 
shock acceleration. The injected number density of relativistic electrons scales 
with the thermal electron energy density, oc pV 2 , rather than the density, p. 
The evolution of the flux is strongly connected to the deceleration of the shock 
wave. The total energy density of the relativistic electrons, if extrapolated to 
7 ~ 1, is ~ 5 x 10~ 4 of the thermal energy density. The free-free absorption 
required is consistent with previous calculations of the circumstellar temperature 
of SN 1993J, T e ~ (2 - 10) x 10 5 K, which failed in explaining the radio 
light curves by pure free-free absorption. Implications for the injection of the 
relativistic electrons, and the relative importance of free-free absorption, Razin 
suppression, and the synchrotron self-absorption effect for other supernovae are 
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also briefly discussed. It is argued that especially the expansion velocity, both 
directly and through the temperature, is important for determining the relative 
importance of the free-free and synchrotron self-absorption. Some guidelines for 
the modeling and interpretation of VLBI observations are also given. 

Subject headings: Supernovae - supernovae: individual (SN 1993J) - radio - 
circumstellar matter - particle acceleration - synchrotron emission 

1. INTRODUCTION 

The radio emission from young supernovae provides an important probe of the 
circumstellar environment of the supernova (Chevalier 1990, Fransson 1994). In addition, 
the time-dependence of the physical conditions makes it a unique laboratory for studies of 
relativistic particle acceleration. 

For the first radio supernovae to be well observed, SN 1979C and SN 1980K (Weiler et 
al. 1986) a standard model to explain their radio light curves was soon developed (Chevalier 
1982b). The basic feature of this model is external free-free absorption in a circumstellar 
medium of the synchrotron emission, arising close to the supernova shock wave. Because of 
the expansion, the optical depth decreases as Tg oc A 2 R~ 3 oc A 2 t -3 , if the shock velocity is 
constant. The intrinsic emission from the shock, i.e. the injection efficiency of relativistic 
particles, as well as the magnetic field, were assumed to scale with the thermal energy 
density behind the shock. This and alternative scalings have been discussed by Chevalier 
(1996). This model has been successful in explaining the observations of 'standard' Type 
IPs (Chevalier 1984), and with a knowledge of the temperature the mass loss rate of the 
progenitor can be determined. Calculations of the temperature (Lundqvist & Fransson 
1988) show that T c « 10 5 - 10 6 K, and mass loss rates of 3 x 1(T 5 and 1.2 x 1(T 4 M yr' 1 
were derived for SN 1980K and SN 1979C, respectively. The model has since been applied 
to a number of other supernovae, although the observational data have for most of these 
been less complete. 

With the recent Type lib SN 1993 J the most complete data set of any supernova 
have become available (Van Dyk et al. 1994; Pooley & Green 1993). When the free-free 
absorption model was applied to this it was, however, found that only a poor fit was 
obtained (Lundqvist 1994; Fransson, Lundqvist, & Chevalier 1996, in the following FLC96; 
Van Dyk et al. 1994). An improved, although still unsatisfactory, fit was obtained if the 
external density instead of p oc r~ 2 followed p oc r~ d where d — 1.5 — 1.7. The reason 
for this flatter profile was not clear, although speculations about variations in the mass 
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loss rate have been proposed. Even more complex models could, not surprisingly, improve 
the fit, but at the expense of a larger number of free parameters without obvious physical 
meaning. A serious problem, pointed out by FLC96, was that the temperature in the 
circumstellar medium was ^ 2 x 10 5 K, resulting in a much smaller free-free optical depth 
than necessary to produce the required absorption. 

In this paper we discuss alternative explanations to the free-free absorption model, 
and show that the radio spectra of SN 1993J are very well explained by a combination of 
synchrotron self-absorption and external free-free absorption, explaining the anomalous 
behavior of SN 1993J, as compared to other radio supernovae. Synchrotron self-absorption 
has been proposed to be responsible for the low frequency cutoff in Type lb supernovae 
(Shklovskii 1985, Slysh 1990), and independently from the present analysis has been 
discussed by Chevalier (1998). Chevalier Sz Dwarkadas (1995) have proposed that the turn 
on of the early radio emission of SN 1987A was determined by synchrotron self-absorption. 
Several of the mechanisms below were discussed by Chevalier (1982b), who, however, for 
good reasons discarded them for SN 1979C and SN 1980K. 

In section |2] we discuss the synchrotron emissivity at a finite density of thermal 
electrons, including the Razin-Tsytovich effect, and in section 3 we apply this analysis to the 
spectra of SN 1993J. We first use simple parameterized models with constant spectral index 
and discuss the resulting fits for models with the three different suppression mechanisms, 
i.e., free-free absorption, synchrotron self-absorption, and the Razin-Tsytovich effect. After 
this we discuss a self-consistent model where the spectral index is determined by the various 
energy loss mechanisms. In section 4 we then discuss the results of this analysis, and 
their implications for the magnetic field, relativistic particle acceleration and circumstellar 
temperature. Section 5 contains models for the full light curves, while in section 6 we 
discuss the relative importance of the Razin-Tsytovich effect, synchrotron self- absorption 
and free-free absorption. In section 7 we make some comments on the analysis of VLBI 
observations of supernovae. The conclusions are summarized in section 8. 



2. ABSORPTION AND SUPPRESSION MECHANISMS FOR THE 

SYNCHROTRON EMISSION 

A main advantage in analyzing the individual radio spectra for the various observed 
dates is that they are only a function of the spectral index of the synchrotron emission, the 
normalization of the spectrum, and a one-parameter or at most two-parameter specified 
function, describing the effect of the external or emitting medium. In contrast, an analysis 
based on the light curves of the different wavelength bands involves uncertain assumptions 



-4- 



about the particle acceleration mechanism (e.g., Chevalier 1996). From the time variation 
of the parameters involved in the fitting of the individual spectra we can instead hope to 
gain some understanding of these uncertain, but extremely important, properties of the 
particle acceleration. 

For a homogeneous shell with a synchrotron optical depth r s (i), and with an external 
free- free medium, parameterized by the optical depth Tg(t) at a wavelength A = 1 cm, 
Gaunt factor gs(X), and impact parameter s, the spectrum is given by 

F„(\) = 2tt S v jf* [1 - e~ T ^) e~ Ts ^ 9ff(A) A " s ds (1) 

where 

"^rp^ CO 

and 

roo Ks (r)n e (r) 2 r 
T * M = h (r* - 8>E*)V* dr - (3) 
Here, Kg is the free-free opacity at 1 cm and n e the thermal electron density. S u is the 
source function, S u = j v /n v , with k u the synchrotron self-absorption opacity and j u the 
emissivity. Finally, R s and D are the radius of the emitting shell and the distance to the 
supernova, respectively. We assume that the supernova ejecta are opaque to the radio 
emission, as is most likely the case. 

In the analysis in section |3.1.1|, we assume that the column density of electrons have 
a constant power-law, giving a synchrotron spectrum with index a. If we further replace 
equation (JJ) by an average over the disk we then get 

F V (X) = S'(t) A" 5 / 2 [1 - e- f ^ A(Q+5/2) ] e" f «« ^ x * (4) 

where the parameter S'(t) is related to S v at 1 cm by S' = tiR 2 s S u (1 cm)/ D 2 . The average 
optical depth over the disk, f s , is given by f s w 2n s AR em , where Ai? em is the thickness of 
the emitting shell. In the optically thin limit we have 



Rs 



2 

2 



F u (\) « 2tt [^j Jv (t)AR era e~ f ^ ^ A . (5) 

In the general case of a plasma with density n e , magnetic field B, and relativistic 
electron distribution N(j) (integrated along the line-of-sight, and assumed to have an 
isotropic pitch angle distribution), the line-of-sight integrated emissivity is given by 
(Tsytovich 1951; Razin 1960; Pacholczyk 1970) 

3u AR em = 1.87 x 10- 23 B sin0 f°° iV( 7 )/(^) F[ J } d 7 (6) 
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where 9 is the pitch angle, v v = 8.98 x 10 3 nl/ 2 Hz is the plasma frequency, and 
v c = i> Bsm9 7 2 is the critical frequency, where u = 4.20 x 10 6 Hz. F(x) is given by 

rod 



F(x)=x / K 5/3 (z)dz (7) 

where K^/^iz) is the modified Bessel function of order- 5/3. The function f(x) = + 
represents the modification due to plasma effects. Similarly, the synchrotron self-absorption 
optical depth is given by 

For a constant power law distribution of electrons, N(^y) = iV rc i 7 _p , one finds 

/ \ cc / i \ poo rp 

j v AR cm = 9.33 x 1(T 24 (1.40 x 1(T 4 A) (B sm9) {a+1) N rcl ^ g(x) F [^T^) x dx 

9 X (9) 

where 

^) = (i + T^r 1/2 , (io) 

Ar x 

and a — (p — l)/2. The critical wavelength, Ar, above which the emission is suppressed is 
given by 

3 -B sin 9 ( n e \ _1 „ . , 

Ar = - = - Q B sin9 cm. (11) 



4 e n e V 1.56 x 10 9 cm -3 

Equation @ contains the pitch angle between the magnetic field and the electron 
velocity. Microscopically, the magnetic field may be random because of turbulence behind 
the shock. In the integrated flux, as measured by the Very Large Array (VLA), one observes 
even a well ordered field from different directions, because of the spherical geometry of the 
shock. We therefore average equation @ over all angles numerically. We also assume that 
the magnetic field and density are constant over the emitting volume. Further, to isolate 
the wavelength dependence from the variation in the normalization, we split the flux into 
one wavelength independent amplitude, A(t), and one spectral function P(X/ Aro,«), 

F„(\) = A(t) P(A/A R0 ,a). (12) 

Here, 

and Aro = Ar(6> = tt/2). The spectral function P(y, a) is given by 

P (y^ a ) = yl r s in#( a+1 ) [°° (i+_^)-V2irU(i + _L)3/2i a .(«-i) da . tZcos0- ( 14 ) 
2 Jo Jo sin 9 x ' sin 9 x 
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P(y, a) is shown in figure |l| for several values of a. The averaged and non-averaged 
spectrum only differ marginally from each other. The averaged spectrum has a slightly 
shallower slope close to the maximum. 

At wavelengths longer than Aro the emission is severely suppressed. This effect has 
been referred to as alternatively the Razin-Tsytovich effect or simply the Razin effect. In 
the following we will for simplicity use the latter terminology. The vacuum case is obtained 
for Aro — ► oo when f — 1. For wavelengths longer than Aro the pitch angle averaged 
spectrum falls off as 

/ \ \ ( a -V2) 3 3 /2 A 

F v oc iV rc i B a+1 A° — (15) 
\Aro / 

(Simon 1969). The cutoff due to the Razin effect is therefore less steep than that due to 
free-free absorption. This is important for the shape of the spectrum at long wavelengths. 

The relevant values of n e and B depend on the particle acceleration site. If this is in the 
shocked circumstellar gas, n e will be a factor of four higher than the electron density in the 
wind. In FLC96 and Fransson (1998) the X-ray emission of SN 1993 J at early and late time 
is modeled. To reproduce the observed X-ray flux a mass loss rate of M » 5 x 10 -5 M Q yr^ 1 
is required, for a wind velocity of u w = 10 km s _1 . With a factor four compression, and 
X(He)/X(H) = 0.3 (e.g., Shigeyama et al. 1994), the electron density immediately behind 
the circumstellar shock is 



n P = 4.36 x 10 8 M 



10 15 



cm 



1.46 x 10 s M ( ^ T ) I \ cm~ 3 (16) 

V2xl0 4 kms-V V 10 days/ V ; 

The similarity solution by Chevalier (1982a) shows that the density behind the circumstellar 
shock is fairly constant (see also figure 3 in FLC96), and we will use equation ( |i~6"D 
throughout the shocked region. We have here for future convenience defined a normalized 
mass loss, A4, rate by 

^ = (5 x 10- 5 M Q yr-^j (lO km s" 1 ) ^ 

In FLC96 it is argued that V ~ 2 x 10 4 km s _1 at ~ 10 days. VLBI observations (Marcaide 
et al. 1995a,b, 1997; Bartel et al. 1994) give an expansion velocity for the radio emitting 
plasma of ~ (1.8 — 2.4) x 10 4 km s _1 during the first ~ 100 days. Consequently, we scale 
our results to V = 2 x 10 4 km s _1 . The motivation for this and the evolution of the velocity 



is discussed further in section 3.2.3 
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With this density Ar is given by 

A,o = 10.7 M (A) ( _ - -I- s i )' (_L-) cm . , 

The magnetic field in the emitting region therefore has to be less than ~ 1 Gauss, or 
alternatively, M/u w large, for the Razin effect to be important. 



3. APPLICATION TO SN 1993J 

Although radio observations are available for a large number of supernovae (e.g., Weiler 
et al. 1998 for a review), SN 1993 J is unique in terms of both frequency of observations, 
small errors in the observed fluxes, and the large wavelength coverage. It is therefore by far 
the object best suited for a detailed analysis. Observations have been published by Van 
Dyk et al. (1994) and Weiler et al. (1998) for 1.33 - 21 cm from 5 - 1300 days, and by 
Pooley & Green (1993) at 2 cm from 8 - 114 days. In addition, there is one set of VLA 
observations at 923 days from 1.3 - 90 cm (Montes et al. 1995). A few mm observations 
with IRAM and OVRI also exist (Radford 1993; Phillips & Kulkarni 1993a, b). The fluxes 
of the mm observations do, however, not agree at similar epochs, and it is difficult to judge 
the accuracy of these observations. For the modeling of the spectra in this paper we use the 
VLA observations by Van Dyk et al. and Weiler et al. , as well as the day 923 observation 
by Montes et al. In section |5] we also include the Ryle telescope observations by Pooley & 
Green (1993), as well as the VLA dates having only a partial spectral coverage, in particular 
before ~ 10 days. 

Our strategy is first to investigate the simplest models, based on a pure power law 
electron and synchrotron spectrum, studying one suppression mechanism at a time. By 
analyzing the results of this, in particular the consistency of these models, we test their 
validity as well as explore the likely range of the relevant physical parameters, in particular 
the magnetic field. However, as a result of the failure of these simple models we are in 



section |3T^ led to a more detailed analysis, taking all the relevant physical processes into 
account. 



3.1. Constant spectral index fits 
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3.1.1. One component models 

We will now use equation (H), assuming here a constant spectral index a = (p — 1)/2, in 
the three limiting cases: 1) free-free absorption, r s = 0, Tg ^ 0, A <C Aro; 2) synchrotron-self 
absorption, r s 7^ 0, = 0, A <C Aro; 3) Razin suppression, r s = 0, = 0, A ~ Aro. In 
the free- free and Razin cases we use the optically thin expression, equation @. In all three 
cases there are three parameters characterizing the spectrum: 1) the overall normalization 
of the synchrotron spectrum, 2) the power law index of the radiation (a), 3) the optical 
depth to free-free absorption (rg), synchrotron self-absorption (r s ), or the Razin cutoff 
wavelength (Aro). For each date we fit the observed spectrum with the set of parameters 
that minimizes \ 2 ■ The value of a is determined by minimizing the scatter of the overall 
light curves; in the Razin case, especially the spectral region where A <C Aro, and for the 
free- free and synchrotron self- absorption cases the spectral region where r(A) < 1. As we 
discuss below there are reasons to expect that, because of cooling, a varies with time and 
frequency. This will be discussed in detail in section [5|. 

In figure |2] we show the observed fluxes at 1.3, 2.0, 3.6, 6.2, and 21 cm at 11.5, 22.3, 
40.1, and 75.1 days. These dates were picked to cover well the exponentially damped part 
of the spectrum. The solid line represents the best fit Razin spectrum, while the dashed 
line is the best fit free-free absorption spectrum and the dotted line shows the same for the 
synchrotron self-absorption case. In the Razin case a = 1.25 ±0.1 gives the best fit. In the 
free-free case a = 0.8 is indicated, while the best fit in the synchrotron self-absorption case 
is obtained for a = 0.7. The reason for this difference in a is that the Razin effect affects 
the spectrum well below Ar , and flattens the overall spectrum. The free-free absorption 
sets in more abruptly, and has less an effect at short wavelengths (cf. equations [| and [TBI) . 

From figure |2] it is obvious that the Razin spectrum gives a superior fit to the 
observations compared to free-free or synchrotron self-absorption. Only when all 
wavelengths are close to the non-suppressed level, as at 40.1 days, do free- free and 
synchrotron self-absorption provide acceptable fits. Otherwise the free-free spectrum 
underestimates the flux, while the synchrotron self-absorption overestimates above the 
wavelength cutoff. At t ^ 200 days all models give reasonable fits because of the decreasing 
importance of both the Razin effect, synchrotron self-absorption, and free-free absorption 
at A ^ 20 cm. 

Because the Razin fit provides the best one component fit, we now analyze the 
implications and consistency of this model in more detail. We have therefore calculated 
the best fit for all dates which have enough frequency points (three for t < 16 days, four 
for 16 < t < 75 days, and five for 75 days < t). The Razin wavelength, Ar , and the 



normalization parameter A(t) (in mJy) (see equation [TJ]) are calculated for these dates 
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and shown in figure |^. During the whole period the Razin wavelength shows a very 
smooth evolution. The solid line in figure |3] shows a simple power law fit to the evolution, 
determined by a least squares fit. Using equation QTTD, the ratio B/n e can be fit with the 
function 

R / + \ a86 

- = 6.6xl(T 10 — ^ (19) 
n e VIOdays; v ; 

with an r.m.s. of only ±0.013 in the exponent. Inserting equation flTE| ) into equation (|T§P, 
the magnetic field can be obtained for the same region as a function of time, 

B^9.6x 10- 2 M ( ^ T ) ( \ ) G. (20) 

V2xl0 4 kms-V \ 10 days/ v ' 

To calculate the total number of relativistic electrons we eliminate in equations (|i~2D 
and fllBD the magnetic field and electron density by using equations (|20|) and (|16D and find 
with p = 3.5 

F, = 1.14 x 10- M 22b ( 2x J kma S * 



X.I t ^°' 5 



F v = 1.6 x 10 2 P(^-) ( \ ) mJy. (22) 



The observed spectral evolution can for t ^ 100 days be fitted by 

\ 0.62 

A R0 ' V. 10 days 
Therefore, 

N ml = 1.4 x 10 25 M- 2 - 25 ( ^ T ) 25 ( —t—) cm' 2 (23) 

V2xl0 4 kms-V VlOdays/ V ; 

for t < 100 days. 

As a consistency check of the assumption of r s <C 1, from equation (§) we can now 
calculate the synchrotron self-absorption optical depth with the above values of B and iV re i. 
One then finds that, due to the Razin suppression of the opacity, the optical depth at Aro 
does not exceed 0.1 at any time. 

A fatal problem for the Razin model occurs when we estimate the energy density 
of relativistic electrons, u Te \. We will be conservative here and just include electrons 
with Lorentz factors responsible for the observed spectrum between 10 - 100 days, i.e., 
wavelengths 1.3 - 6.3 cm, giving 1.1 x 10 2 (t/10 days) - 57 < 7 < 2.4 x 10 2 (i/10 days) - 57 . 
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From equation fl23|) and p = 3.5 we find that the energy per unit area within this band 
is ~ 4.7 x 10 15 (t/10 days) -0 ' 75 ergs cm -2 , which is only weakly dependent on p. As 
we show in section |3.2.1| (equation [TF] with L bol rj 4 x 10 42 (t/10 days) -0 - 9 erg s -1 ), 
inverse Compton losses are in this case likely to be important down to 7 ~ 30. This 
implies in the Razin case that the inverse Compton timescale during the first 100 
days is t Comp ~ 0.36 (t/10 days) 2 ' 4 A 0,5 days. The thickness of the emitting region is 
AR em « Wcomp/4, or for A ~ 3 cm Ai? cm w 2.5 x 10 13 (t/10 days) 2 4 cm. The energy 
density is then u re \ ~ 1.9 x 10 2 (t/10 days) -3 ' 1 ergs cm -3 . This can be compared to the 
thermal energy density behind the shock, ^thermal ~ 3.4 x 10 2 M. (t/10 days) -2 ergs cm -3 , 
which is only slightly higher than u Te \. 

Our estimate of u re \ is most likely a serious underestimate, since we have only included 
the spectrum between wavelengths 1.3 - 6.3 cm (already at 21 cm, 7 ~ 60 (t/10 days) ). 
In reality one should extrapolate the electron spectrum down to 7 ~ 1, which could increase 
the energy density by a factor of j(X = 1.3 cm) 1 ' 5 ^ 300. The exact factor depends on 
the detailed cooling processes as well as the form of the spectrum close to 7 ~ 1, and is 
therefore model dependent. This will, however, not alter our main conclusion that the 
energy density of relativistic electrons in this model is likely to exceed the energy density 
of thermal particles by a large factor, which is unrealistic in any model for the particle 
acceleration, and rules out this model. Typical values of the injection efficiency are instead 

10 -2 (e.g., Chevalier 1982b). The main reason for this problem is the weak magnetic 
field and the steep power law index derived from the Razin cutoff. Because F v oc iV rol B a+1 
this results in a very large value of N Te \. 



3.1.2. Synchrotron self- absorption plus free-free external absorption 

The fact that the Razin model requires an unrealistically large density of non-thermal 
electrons, due to the low magnetic field, forces us to more complex models. The most 
natural one is then a combination of synchrotron self-absorption and external free-free 
absorption. The latter component has the effect of suppressing the too high flux at long 
wavelengths in the synchrotron self-absorption model. As we will find, this model also 
solves the efficiency problem above by requiring a higher magnetic field, while at the same 
time giving a free-free absorption consistent with earlier estimates in FLC96. 

To model the combined absorption we now use equation (|l|), integrated over the 
supernova disk, but still assuming a fixed spectral index. We then have r s (l cm) and 
Tff (1 cm) as free parameters, in addition to a and S'. 
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For a given temperature, T e (r), and including the Gaunt factor, g ff , the free-free optical 
depth along a ray with impact parameter s is given by 



,53 ,',2,2 U + 2A') r (-2.78 + log[A TV 2 (r)]) 



Mt, s) = 5.63 x 10" AW ^_ 7 _HL_^ ^ (24 ) 

where X = X (He) / X (H.) . Tg is therefore sensitive to the radial temperature profile, as well 
as its dependence on time. Lundqvist & Fransson (1988), as well as FLC96, have discussed 
the determination of this extensively, and it is found that at radii, ^ 1 x 10 16 cm the gas 
cools rapidly from > 1 x 10 6 K to ~ 2.5 x 10 5 K, while at large radii the temperature is 
fairly constant at ~ (2 — 2.5) x 10 5 K. Therefore, as the shock propagates outwards, the gas 
close to the shock will initially have a decreasing temperature, while later it will be more 
constant. This is seen clearly for SN 1993J in figure 11 of FLC96. While this behavior is 
probably quite generic, the details of the radial variation depend strongly on the properties 
of the shock radiation (see section 14.11). 



In this paper we use a simple parameterization of the temperature profile, assumed to 
be constant with time, and given by 



Tp(r) = max 



(25) 



We then vary the parameters T15 and 5 to get an optimum fit of the spectra. Although 
indicative of the actual temperature, equation should only be interpreted as a 



convenient fitting formula. The reason is that the fundamental parameter for the fitting 
of the spectra is not T e (r), but rather Tg(t). As shown by equation (|24"|), the same rg(t) 



can be obtained by different Ai and T e (r), as long as Ai 2 J ' T e (r)~ 3 / 2 r~ 4: dr » constant. 
Different temperature profiles can therefore yield the same result for Tg(t). The value of T e 
given by equation (|25|) should therefore only be interpreted as a weighted average. Also, 
the dependence of T e on Ai should be kept in mind when interpreting the results. 

We have calculated the best fit for different values of a, T15, and 6, varying S'(t) and 
r s (t) for each date. We find that good fits are obtained for 5 = 1.0 and T 15 = 2.75 x 10 6 
K, with xt = 1-4 — 1-5. Larger and smaller values of 5 give considerably higher xt- The 
free-free absorption is only well-determined by the spectra and light curves later than ~ 10 
days. We have tested different values of the spectral index, and find that we get a best fit 
with a = 0.85, or p = 2.7. Because the quality of the fit is very similar to that in the more 
detailed model in figure |8], we do not show it here. 

In figure |] we show the evolution of the wavelength at which the synchrotron 
self- absorption optical depth is unity, A Tssa= i = r s (t) _1 /( 5//2+a \ and the scaled source function 
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at 1 cm, S'(t). A least squares fit to the observations for 10 < t < 100 days gives 

with Aio = 1-07 cm, and q\ = 0.68. For t }t 100 days the evolution gets steeper and 
Aio = 0.81 cm, and q\ = 0.81. The dispersion in the exponent is ±0.03. Similarly, the 
variation in S' can for 10 - 100 days be fitted by 



10 days 



S'(t) = S 10 ( ] , (27) 



where Sw = 60.1 mJy and qs = 2.43, while at later times ( }t 100 days) the evolution is less 
steep with Sw = 1.67 x 10 2 mJy and qs = 2.00. 

Knowing A Tssa= i(t) and S'(t), we can now calculate the magnetic field and the column 
density of relativistic particles, 

/ p \ 4 
,2 / J^s \ o'/j-\-2 



Bit) = 9.45 x 10- 8 c(p) 2 (^J S\t)~ 2 , (28) 
where 

r(n) _ (p+j) r(g±g)r(^)r(^)r(3gi) 

[P) (p + f)(P + 1) r(£?)r(2±Z)r(^)r(^) [ZJ) 

and T is the gamma function. The column density of relativistic particles is then given by 

N (+\ 1.98 X 10 17 rl 2 ^) n-( P +2)/2 r (f+4)/2 „ m -2 (nn\ 

(i.i8xio- 2 )Mp+f)r(^)r(^)r(3£±iO) T - =1 ' 1 } 

As is seen from equations (^) and (|3!1D, both B{t) and N re i(t) depend on R s (t), which 
in turn is determined by the ejecta structure. In this section we assume for simplicity a 
constant velocity of V = 2.2 x 10 4 km s -1 (see section |3.2.3| ). 

The resulting magnetic field and column density are given in figure [5[ During most of 
the observed interval the decrease in B is approximately a power law, and such a fit gives 

/ ^ \ -0.93±0.08 y x 4 



Bit) = 25.5 G. (31) 

w VlOdays/ V2.2 x 10 4 km s" 1 / v ; 

Within the errors, N ie i(t) is nearly constant, although the scatter increases considerably 
for t ?t 100 days. The level of iV re i in figure [| should be compared with that in the Razin 



model from section |3.1.1|. Because of the higher magnetic field, and most importantly the 



smaller value of p, N re \ is roughly six orders of magnitude smaller than for the Razin model. 
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With the magnetic field in equation (|31| ) we can calculate the Razin wavelength from 
equation (|T8| ) 



A R0 = 3.5 x 10 2 M ( \ r ) ( - ] 

V2.2 x 10 4 km s- 1 / \10 days / 



Comparing equations (|26|) and (^) we find that the Razin effect in this model is 
unimportant. 

Although the above results are of great interest, and will be discussed in detail below, 
there are several deficiencies in the analysis. The most serious limitation is that we have 
assumed a constant power law electron spectrum, i.e., neglected possible losses by radiation 
or collisions. We will therefore now, at the expense of a more complicated model, discuss 
a self-consistent analysis of the observations. The discussion in this section is, however, of 
great help in order to understand the different processes involved, as well as a guide to the 
values of B and N Te \ and their sensitivity to different parameters. 



For formulating a self-consistent analysis it is necessary to discuss possible energy 
losses of the electrons which will affect the integrated electron spectrum, and therefore also 
the emission from the plasma. To estimate these effects we assume in this subsection that 
the magnetic field behaves like 



where B w parameterizes the strength of the field. From our earlier discussion and 
anticipating the results below, we expect B w pa 30 G, to which we scale our estimates. 
To express our results as a function of time rather than radius we use in this section 
R(t) = 1.72 x 10 15 (t/10 days) cm, corresponding to a constant velocity of V = 2 x 10 4 km s -1 . 

We first estimate the typical Lorentz factor for the relativistic electrons responsible for 
the radiation at wavelength A from 7 pa 85 (A B)^ 1 ^ 2 , or 



3.2. Self-consistent analysis 



3.2.1. Energy losses of the non-thermal electrons. 




(33) 




(34) 



The synchrotron lifetime for an electron to decrease its energy by a factor two is 




(35) 
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In terms of the observed wavelength this becomes 



^ = 6.4 x lO" 2 A 1/2 ( ) — i— (36) 



t V30 G/ \ 10 days, 

It is obvious that synchrotron losses are important throughout most of the observed period. 

Because of the strong background radiation field from the supernova, inverse Compton 
losses may be important for the relativistic electrons. One can estimate the timescale for 
this from 

W = A «4q 7 -i ( L Boi rV V \ 2 ( _t_ ) (37) 

t 167rja T J Bol 1 \10 42 ergs- 1 / V2 x 10 4 km s -1 / \^ 10 days/' v ; 

where Jb i and L Bo i are the bolometric mean intensity and luminosity, respectively. 

The bolometric luminosity has contributions from both the ejecta, powered by the 
thermal energy from the passage of the shock and later by radioactive decay, as well as 
radiation produced by the circumstellar interaction. The ejecta flux is emitted mainly in 
the optical, and at early time also in the UV. We estimate the ejecta luminosity from the 
bolometric luminosity calculated by Young, Baron & Branch (1995), from observations by 
Richmond et al. (1994) and Lewis et al. (1994). Unfortunately, the reddening is uncertain. 
A range in E-q-v from 0.08 - 0.32 has been estimated by Richmond et al. Unless otherwise 
stated we will use Eb-v = 0.08 in the following. 

The luminosity in X-rays from the circumstellar and reverse shocks are calculated 
in FLC96 and Fransson (1998). The circumstellar shock is important only for 
t ^ 10 days, and based on figure 8 in FLC96, scaled to M. = 1, we estimate 
L cs = 3.1 x 10 41 (t/10 days) -1 ' 2 erg s _1 . To include both the in-going and outgoing radiation 
we multiply the results in FLC96 by a factor two. In Fransson (1998) the evolution of 
the reverse shock luminosity is calculated for realistic density profiles, and it is found that 
L rev = (1 - 2) x 10 41 erg s _1 . 

The mean intensity, J,, of the different luminosity components at the position of the 
circumstellar shock is given by J 4 = WLi/A^R 2 , where W — |[1 — (1 — (i?j/i? s ) 2 ) 1//2 ] is the 
dilution factor. The radius Ri stands for either the photospheric radius or the shock radius 
in the two cases. For R p h -R s , J ~ L^/1Qtt 2 R 2 . The radiation from the shocks is more 
isotropic at the circumstellar shock. Depending on the thickness of the circumstellar shock 
and the location of the acceleration zone, one expects that R s ~ (1.0 — 1.3) -R rev - 

In figure ^ we show the resulting total mean intensity multiplied by 167r 2 i? 2 , together 
with the individual contributions from the ejecta and the shock waves for M. = 1. In the case 
that Ri -C -R s this essentially gives the total luminosity from component i. During the first 
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~ 100 days the radiation density is dominated by the ejecta, with an increasing contribution 
from the reverse shock, and can be approximated by L^ Q \ « 4 x 10 42 (t/10 days) -0,9 erg s -1 . 

The Compton scattered optical photons from the ejecta and the thermalized X-rays 
from the reverse shock emerge at hv « 7 2 hvo ~ 0.1 — 100 keV. For our parameters the 
luminosity in this range is ~ 2.5 x 10 35 (L synch /10 37 erg s _1 )(L bo i/10 42 erg s" 1 ) erg s _1 . The 
free-free flux from the shock (FLC96) therefore dominates the X-rays. One can also show 
that the ratio between first- and second-order Compton scatterings (Rees 1967) is small. 

Knowing both the magnetic field and the radiation density, we can calculate the critical 
value of the Lorentz factors, 7 syn ch and 7com P , for which the synchrotron and Compton 
timescales are equal to the expansion timescale. In figure |7| this is shown as the solid line. 
We also show the approximate value of 7 for the electrons responsible for the emission 
at 1 cm and 21 cm, 71 cm and 721 cm , respectively, as given by equation (|4]). Clearly, 
synchrotron losses are more important than inverse Compton losses at all times. 

Coulomb losses can also be important for the relativistic electrons. For these electrons 
the ratio of the Coulomb timescale to the expansion timescale is given by 

^ = 8.9 x 10- 3 7 M- 1 ( £ =V ( —t—) • (38) 

t 1 V2xl0 4 kms-V \ 10 days/ v ; 

For 7 ^ 7coui = 11 x 10 2 (t/10 days) -1 Coulomb losses are important. We should also 
compare tcoui with the synchrotron timescale, and we find that for 

7 < % -^8.62^ (||)- 1 ( -_^_ r )- 1 (39) 

Coulomb losses dominate synchrotron losses. In figure [7|we show 7c ou i an d 7 s -c for M = 1. 
We note that 7 s _c is independent of time, and that especially at early times Coulomb 
cooling is important for the longer wavelengths. 

Synchrotron and inverse Compton losses steepen the electron spectrum index by one 
(e.g., Pacholczyk 1970, and section [| below). Coulomb losses have the opposite effect to 
synchrotron losses in that the electron spectrum index flattens by one power, compared to 
the injected spectrum. Therefore, with the parameters in figure ff] we expect that before 
~ 100 days (i.e., when 7 sync h ~ 7coui) the electron spectrum will be flattened by one power 
for energies below 7 s _c and steepened by one power above 7 s -c- In reality, both Coulomb 
and synchrotron compete in the regions responsible for the observed emission, and the 
spectral index is intermediate between these values. Later than ~ 100 days the electron 
spectrum will have three regions, one between 7c ou i < 7 < 7s y nch, where the spectrum is 
equal to the original input spectrum, while below and above this region the spectrum is 
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flattened and steepened by one power in energy, respectively. We return to this issue in 
section ||. 



3.2.2. Calculation of the electron spectrum 

For the acceleration of the relativistic electrons, once they have been injected from the 
thermal pool, there are at least two possibilities. First-order Fermi acceleration across the 
circumstellar shock is an obvious candidate, but acceleration in the turbulent region close 
to the contact discontinuity is also a possibility (Chevalier, Blondin, & Emmering 1992). 
The timescale of acceleration depends on the spatial diffusion coefficient, and therefore 
on the mean free paths. Only rough estimates of the timescale based on modeling of the 
radio emission from SN 1987A are available (Ball & Kirk 1992; Duffy, Ball, & Kirk 1995). 
From the expression given by Ball & Kirk (1992) for the variation of flux with time (their 
equation 8) we find that the rise is likely to be very rapid, t acc <C 10 days. Here, we therefore 
assume that the electrons are accelerated almost instantaneously (i.e. t acc <C tcomp,^), with 
a spectrum cx 7 _Pi . 

In calculating the electron distribution we use the equation 

t - fa™ + - ^ + ^ p?r • w 

The first term on the right hand side takes all energy losses into account, with a total 
energy loss timescale given by 

/ 1 1 1 IV 1 , . 

tloss = + + + - • (41) 

\ ^synch ^Comp ^Coul « / 

The first three terms on the right hand side are already defined, while the last term takes 
adiabatic losses into account. In equation (|40|) the second term, where Efa) is given by 

£ w=4?/ j(!/ ' 7) ^ (42) 

accounts for the heating of the non-thermal electrons by the radiation. In this 7) is 
the single electron emissivity. The third term describes the escape of particles, either as a 
result of advection towards the contact discontinuity, and subsequent energy degradation by 
collisions in the high density gas, or as a result of actual escape. Finally, the last term gives 
the number of injected electrons with energy 7m e c 2 per second and unit area, which we 
assume to be instantaneously accelerated at the shock. A similar equation was first used by 
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McCray (1969) to discuss the thermalization of the electron distribution, when synchrotron 
self-absorption is important, while the form given here is basically that used by de Kool, 
Begelman, & Sikora (1989). The equation assumes that the electrons are relativistic, but 
given our ignorance of the details of the injection process, which certainly influences the 
distribution for 7 ~ 1, this is sufficient. Equation (|40f) is solved by the semi- implicit method 
of Chang & Cooper (1970). 

If we assume that the self-absorption effects are small and that the escape term can 
be ignored to a reasonable approximation, as is confirmed by the detailed calculations in 
section |5|, we can write down the stationary solution to equation ( |40| ) as 



tf(7) = -£— -M ( 1 + 7— + 7 + 7 ) ( 43 ) 

^ ;min \ /min / ^Coul ^Comp ''synch 

With this electron spectrum we can calculate the synchrotron emissivity and opacity from 
the general expressions in equations © and (§), as well as the source function, and finally 
the synchrotron spectrum as function of time. 

Because V(t), Lboi and M are determined from other observations (section |3.2.3| ), we 
have for each spectrum two parameters, n re i and B. We have to assume a specific input 
spectrum for the electrons (i.e., pi), which we take to be constant with time. 



3.2.3. Velocity evolution 

The magnetic field and electron density depend on the velocity evolution of the shock, 
i.e., on R s (e.g., equation 28). The velocity of the shock region is in turn sensitive to 



the density structure of the ejecta. If the ejecta density is approximated by a power law 
with an exponent n, the shock radius varies as R s cx t m , where, for an n e oc r~ 2 wind, 
m = (n — 3)/(n — 2) (Chevalier 1982a). A flat density profile, n < 7, causes the shock 
to slow down, while a steep profile, n ^ 20, gives a shock of fairly constant velocity. 
Unfortunately, the density structure of the ejecta is uncertain, depending on an equally 
uncertain progenitor structure. Optical line profile observations of SN 1993J provide some 
constraint on V(t). This, however, refers to the ejecta gas, and not directly to the radio 
emitting plasma, which is likely to be 20 — 30 % further out. Also, if there is a change in 
the density gradient, one can expect a lag of the adjustment to the new gradient. 

The most direct information about the size of the radio emitting region comes from 
VLBI observations. During the first 92 days Bartel et al. (1994) found roughly undecelerated 
expansion, m = 0.96 ± 0.07, with a velocity V = 19000 ± 3000 km s -1 . Most of the 
uncertainty in V comes from the M81 distance errors, D = 3.63 ± 0.34 Mpc (Freedman et 



-18- 



al. 1994), which do, however, not enter in m. At later times, 182 — 1304 days, Marcaide et 
al. (1997) find strong evidence for deceleration, with m ~ 0.86 ± 0.02. It is quite possible 
that the VLBI velocities above are underestimated. FLC96, as well as Bartel et al. , argue 
that the Ha line indicated an ejecta velocity of 18000 - 19000 km s _1 at 10 — 30 days. The 
shock velocity should be a factor of 1.2 — 1.3 larger than this, or 21600 — 24700 km s _1 . 
Furthermore, the VLBI measurements have a fairly small dynamic range; faint high velocity 
material may therefore go undetected. We also point out that the radius measured by 
Marcaide et al. only refers to the 50 % contour level, which explicitly excludes such faint, 
high velocity material. Guided by this discussion we will in the following parameterize the 
velocity by 

(V t<t p 

with Vq = 22000 km s _1 . Because we take V = constant for t < t p our value of m will be 
less than the value of m = 0.86 from Marcaide et al. With t p = 100 days our effective value 
of m up to day 1304 is m = 0.74, corresponding to n fh 6. Using t p = 200 days would not 
change our results significantly. 



4. RESULTS 

With all parameters fixed, except B, n re \ and p- u we can now repeat the modeling in 
section |3.1.2| for each date, thereby determining the unknown parameters from a x 2 fit- 



However, in contrast to section |3.1.2j we take the different cooling processes into account. 
The resulting spectra at six epochs, from 11.5 - 923 days, are shown in figure |8|. Also 
shown is the spectrum when free-free absorption is neglected. It is seen that when the latter 
process is included the fit is essentially perfect, as is the case also for the other dates. We 
now discuss the implications for the physical parameters resulting from this procedure. 



4.1. Free-free absorption 



The analysis in section |3.1.2| shows that the best fits for the free-free absorption are 
obtained for 5 = 1.0 and T15 = 2.75 x 10 6 K and an asymptotic temperature of 2 x 10 5 K. 

This, and the evolution of temperature derived from the observations, can be compared 
to the calculations in FLC96. At 10 days FLC96 (their figure 11) find that at a radius of 
1.5 x R s , T e m 6.3 x 10 5 K, decreasing to 3.0 x 10 5 K at 20 days and then staying at a 
constant value of 2.5 x 10 5 K up to day 100. 
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As we stressed in section 3.1.2 , the most relevant quantity to compare with the 
observations is not the temperature, but the the optical depth. A further complication 
is that the temperature in FLC96 not only varies with radius, but, contrary to our 
model assumptions, it also shows some variation with time at a given radius. For a real 
comparison one should perform the radial integration of the optical depth at each time (i.e., 
r ff oc jR s T e {r,t)~ l ^n e {r) 2 dr). In figure |9]we show the free-free optical depth, averaged over 
the disk, for the p oc r~ 2 model in FLC96 along with that determined by the observations. 
In general the two results agree to better than a factor of two. While the difference is 
significant, we feel that it is remarkable that they agree as well as they do, and that 
this agreement gives us increased confidence in the determination of the circumstellar 
temperature from the model calculations not only in the case of SN 1993J, but also for 
other objects like SN 1979C and SN 1980K (Lundqvist & Fransson 1988). An approximate 
fit to our results in figure is given by 

*w = ai2A "(id^)""- (44) 

The fact that we get a power law significantly flatter than t~ 3 even in the constant 
temperature region at t ^ 100 days is due to the compensating effect of the decreasing 
velocity, which gives an approximate power law dependence over the whole period. 

Although in broad agreement with the optical depth required by our modeling, there 
are considerable uncertainties in the calculation of the temperature of the circumstellar 
medium, as is discussed in FLC96. Most important of these are the degree of equipartition 
of electrons and ions behind the circumstellar shock, and the evolution of the shock velocity 
shortly after break-out. Both affect the temperature of the shocked electrons, as well as 
the electron scattering optical depth. These two quantities determine the UV and soft 
X-ray flux from the shock through Comptonization of the soft photospheric photons. The 
uncertainty in the shocked electron temperature is therefore reflected in the absolute value 
and the radial dependence of the temperature of the circumstellar medium. The radio 
observations therefore provide a new constraint on the physics and radiation mechanisms of 
the shock. 



4.2. Magnetic field 

In figure we show the resulting magnetic field as a function of the shock radius. 
The errors in B and n re \ are determined from the 1 a contours of constant \ 2 m the B - 
n re \ plane. The scatter in B{t) at late times is related to the scatter in S'(t), which enters 
quadratically in B{t). 



-20- 



A least squares fit to the data shows that 



with B 1S = 61.1 ± 4.7 G and (3 = 0.984 ± 0.028. As seen by the dotted line in figure [TOj, it is 
certainly compatible with a simple B oc R~ l law. A best fit with this dependence is given 
by .Bis = 63.5 ± 4.8 G (f3 = 1.0). Note that the errors do not include systematic errors. 
In particular, the results are, as shown by the expressions above, sensitive to the velocity 
evolution, as well as its absolute value. 

A B oc R~ l dependence is close to that expected if the magnetic field behind the 
shock is amplified by a constant factor from the circumstellar magnetic field of a rotating 
progenitor. If the wind has sufficient angular momentum, the magnetic field should have 
the topology of a 'Parker spiral' at large radii (Parker 1958; Weber & Davis 1967). In this 
case one expects B oc r _1 , while a radially expanding wind has B oc r~ 2 . 

The magnetic fields of the circumstellar media of late type supergiants are uncertain. 
Based on polarization observations of OH masers in supergiants, Cohen et al. (1987) and 
Nedoluha & Bowers (1992) estimate that at a radius of ~ 10 16 cm the magnetic field is 
~ 1 — 2 mG, though the uncertainty in this number is large. It is unlikely that the magnetic 
field in the wind is higher than that corresponding to equipartition between the magnetic 
field and the kinetic energy of the wind. This means that ^- % p giving 



10~ 5 M yr- 1 ) VlOkms" 1 / VlO 16 cm 

Likely locations for the electron acceleration are at the position of the circumstellar shock, 
or, alternatively, close to the contact discontinuity between the circumstellar swept-up 
gas and the shocked ejecta gas. The latter region is Rayleigh- Taylor unstable, and the 
associated turbulence may help amplifying the magnetic field (Chevalier, Blondin & 
Emmering 1992; Jun & Norman 1996). 

At 10 days, corresponding to a radius ~ 1.9 x 10 15 cm, we find that the magnetic field 
in the emitting region is ~ 34 G and the . Using the above estimate of the circumstellar 
magnetic field and a shock compression by a factor four, this post-shock magnetic field 
would be B fa (2.4 — 4.8) x 10~ 2 G. This is a factor ~ 10 3 less than that inferred from 
the observations, and therefore strongly argues for magnetic field amplification behind the 
shock. Although this conclusion rests on the very uncertain estimate of the circumstellar 
magnetic fields of the progenitor system, a simple shock compression of the field can 
probably be excluded, and gives support to the results found in the simulations by Jun & 
Norman (1996). 
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If the scaling of the magnetic field is based on the thermal energy density one expects 
B 2 /8n oc pV 2 , implying B oc V/R s oc t^ 1 oc R~ l l m (e.g., Chevalier 1996). This is shown as 
the dashed line in figure [II]. Although there is formally better agreement with the B oc R~ x 
curve (dotted line), the errors in the data in figure |H] are large enough that we find it 
difficult to separate these two cases. In the following discussion we will assume that the 
magnetic field follows 

m = 1 u.»> toi - )' • (47> 

with B w = 34 G. This choice is of no importance for the subsequent analysis, and is only 
argued from the point of view of the possible importance of equipartition. 

The ratio of the magnetic energy density to the thermal energy density behind the 
shock is 



u b n , , ( Bio 



- 0.14 (-^f M- 1 (48) 

independent of time. It is interesting that the field is of the same order as that given by 
equipartition. The difference from equipartition is, however, significant. 



4.3. Injection rate of non-thermal electrons 

As we discuss in more detail in the next section, the power law index of the injected 
electron spectrum is close to pi = 2.1, and constant in time. We therefore use this value in 



the discussion below. In figure |TT] we show the density of the injected non-thermal electrons, 
given by n Te \, as a function of shock radius. The value of n ve \ is determined mainly by the 
optically thin flux, and therefore, from the arguments in section |5|, can be shown to depend 
on velocity as n re \ oc V 3 ~ 2pi oc V^ 1 ' 2 . A least squares fit to the data in figure [IT] for the first 
100 days is given by 

n rel = n rel 15 7m - ( T ^-) ' ( 2 2x J ^^ ) , (49) 

where n re \ 15 = (6.1 ±0.7) x 10 4 cm -3 and rj = 1.98±0.04. There is after 100 days a prominent 
steepening of the slope, and one finds that n re \ i 5 = (4.0±0.9) x 10 5 cm -3 and r\ = 2.64±0.05. 
A fit based on r\ = 2 during the first 100 days gives n ve \ 15 = (6.4 ± 0.8) x 10 4 cm -3 . 

Chevalier (1996) has discussed scalings for the number density of relativistic particles 
based on either a constant fraction of the thermal particle density, n re \ oc p w i n d oc R~ 2 , or 
a constant fraction of the thermal energy density, n re \ oc p w i n d V 2 oc R~ 2 V 2 oc t~ 2 . Here 
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Pwmd is the wind density. These scalings have little physical motivation, and can only be 
justified by observations. An argument in favor of these scalings is that similar scalings 
have reproduced the light curves of other supernovae dominated by free-free absorption. 
We can now test these scalings directly against the observations in the same way as the 
magnetic field. Taking the two alternatives for the relativistic particle injection efficiencies 
into account, we write 

fV\ 2e 

n re \ = /rei n e y—j , (50) 

where e = for n rc \ oc p w i n d and e = 1 for n re \ oc p w i n d V 2 . The thermal electron density, n e , 
behind the shock is given by equation flIE|) . 

As long as the velocity is constant both prescriptions give the same result, n re i oc R~ 2 . 
When the velocity decreases after 100 days, p win d V 2 oc t~ 2 oc R~ 2 ^ m oc R~ 2m for m = 0.74. 
This is shown as the dashed line in figure [11]. As is obvious from both the figure and the fit 
above, n rel oc p W md V 2 reproduces the observations remarkably well, while the n rel oc p w i n d 
model is incompatible with the observed evolution. This important result shows that a fixed 
fraction, f re \, of the thermal particle energy density is accelerated to relativistic energies. 

The best fit value of / re i found by this analysis and the modeling in next section is 
/ r ei = 1.85 x 10 -4 . The energy density in relativistic electrons immediately behind the shock 
is 

(7mfn +2 - 7max +2 ) 

"rei = fre\n e m e c 2 min - ™ x — = frein e m e c 2 Q (j min , 7 max , p ; ) . (51) 

The ratio of the non-thermal electron density and the thermal total energy density behind 
the shock is M re i/^therm ~ 5.9 x 10~ 5 Q. For p x = 2.1, 7 min = 1 and 7 max = oo, Q = 8.7, 
while 7 min = 5 and 7 max = 15, corresponding to the directly inferred range at 10 days, gives 
Q = 0.9. In either case the non-thermal electron density is much smaller than the thermal 
electron density, in contrast to the Razin case. There may, however, be a substantial energy 
density of non-thermal ions accelerated by the shock. For an equal number density of ions 
one then gets Wion/wthcrm ~ 0.11 Q. It is therefore conceivable that, to within an order of 
magnitude, rough equipartition between the non-thermal particles, the magnetic field and 
the thermal energy is present behind the shock. This conclusion rests, however, totally on 
the ad hoc assumption of equal efficiency for accelerating electrons and ions. 



5. RADIO LIGHT CURVES FOR SN 1993J 



The most direct way to compare the model with the observations is through the radio 
light curves at the different wavelengths. By calculating the time evolution we can also 
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include other effects that were not possible for the analysis above. It therefore serves as a 
check on the assumptions of this analysis. In contrast to the analysis in sections fO and 



we now make an assumption about the the behavior of the magnetic field and injection 
efficiency For the injection rate we use equation (pOf), with e = or e = 1. The high ratio of 



WB/^therm (equation HD may be taken as an argument for a scaling ub oc Wtherm, or £> oc 
rather than B oc R' 1 , and the calculations below therefore assume B oc £ _1 . As we show, 
the results are essentially independent of this assumption. For the radiative transfer we use 
equations (JT|) , (H), and (||). In contrast to section [|, we now use the full form of equation 
(|4Ci|), including self-absorption effects on the electron distribution, as well as advection and 
the explicit time-dependence. 

The evolution of the velocity of the shock (through B, n TC \ and R) is important for the 
number of relativistic electrons, and we use the parameterization in section |3.2.3| . In the 



analysis we include all VLA observations, as well as the 2 cm observations by Pooley & 
Green (1993). In particular, this includes several observations from 5-10 days, not used in 
the analyses in sections |3.1.2| and |j. 



The resulting model light curves are shown in figure [12], together with all observations. 
As can be seen, the model reproduces the light curves at the different wavelengths very 
well. The reproduction of the observed slope at early times is a result of 

F„(r ssa > 1) oc R 2 S u e~ Ta oc R 2 B~ 1/2 e~ TB oc t 2 - 5 e~ TS , (52) 

(e.g., Chevalier 1982b). This is considerably shallower than the pure free-free case discussed 
by FLC96. Because the magnetic field is mainly determined by the optically thick part 
of the spectrum, equation fl5"2"D illustrates the sensitivity of the derived value of B to the 
free-free optical depth. 



In figure [13] we show the electron distribution at three different epochs, and in figure 
|14| the effective spectral index of the electron distribution, p e s = 2a + 1, at the energies 
responsible for the emission at the different observed wavelengths. At early times the range 
is very large, p e R = 1.3 — 2.5, while at t ^ 1000 days the range is smaller at p e R = 2.4 — 2.8 
(yet still higher than p\) , showing that cooling is important throughout the whole period. 
The range of p e g even at late times clearly illustrates the necessity to do a self-consistent 
calculation of the spectrum, as has been done here, rather than an analysis based on a 



constant spectral index as in section |3.1.2| . The evolution of the electron spectrum can be 
compared directly to figure [j] where it is seen that for the shortest wavelengths synchrotron 
cooling dominates, while at longer wavelengths Coulomb cooling dominates before ~ 100 
days, explaining the low value of p e g especially for 21 cm. It is, however, important to note 
that the effects of the various cooling processes extend over a considerable energy range, 
and the values of 7 sync h and 7c ou i are only rough indicators of the cooling effects. 
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Synchrotron self-absorption effects on the electron distribution are modest and can be 
seen only close to the peak in the spectrum (figure [13]). This is the cause of the crossing of 
the spectral index curves seen in figure |TJ]. The reason for the modest self-absorption effects 
on the electrons is that the Coulomb timescale is shorter than the synchrotron timescale for 
7 corresponding to t(v) ^ 1. The electrons are therefore down-scattered in energy faster 
than they are redistributed by the self-absorption, in agreement with what de Kool et 
al. (1989) find. The optically thick flux is quite independent of assumptions about spectral 
index or cooling, which are important when r s ^ 1. Only close to the peak in the spectrum, 
where r s ~ 1, can some effects of the self- absorption on the electron distribution be seen. 

Free-free optical depth unity is reached on days 12.6, 26, and 80 for 3.1, 6.2 and 21 
cm, respectively. At the shorter wavelengths it never exceeds unity because of the high 
circumstellar temperature. The Razin effect is unimportant for the flux at all wavelengths 



because of the necessarily high magnetic field (equation |32]) 



Both the optically thin luminosity and the optical depth are determined by the 
line-of-sight integrated number of relativistic electrons per energy, N(j), and a qualitative 
discussion of this quantity is of interest. For this purpose we neglect the self- absorption 
effects on the electron spectrum, as well as the Razin suppression. If we write equation 
in terms of the synchrotron timescale t sy nch = 67rm e c/<7T7-B 2 , we get 



iV( 7 ) 




t 



10 days 



where 



C{t) = 2.49 x 10" 



B 



10 



30 G 



V 



(53) 



(54) 



^lO 42 erg s- 1 / V2 x 10 4 km s" 1 
The term C(t) gives the ratio of inverse Compton to synchrotron losses, and can in our 
case be neglected. The third term in the bracket gives the Coulomb to synchrotron losses 



(equation [39]), and the last, the adiabatic losses. When synchrotron losses dominate the 
factor 7~(P i+1 ) means that the integrated electron spectrum is steepened by one power, 
while if Coulomb losses dominate it is flattened by one. 



Using equations (|53|), divided by 47r, and a delta function approximation 
for the total emissivity of a single electron in an isotropic magnetic field, 
j(7, v) = c(TtB 2 ^ 2 5{v — z/ i?7 2 )/67r and j u AR era = J N^/An z/)c?7, in equation 
(H) with f ff 1, we find 
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(56) 
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With equations (|T 
i^(r ssa <C 1) 



and (^) we get 
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The slow increase in the flux of the optically thin light curves at 1.3 and 2 cm before 100 
days is caused by the decreasing importance of Coulomb losses with time. Furthermore, 
taking a = d In F v j dlnu and using p e g = 2a + 1 gives a good representation of the evolution 



of the electron spectral index at the various wavelengths, as shown in figure [14 



After ~ 100 days the shock velocity decreases according to m ~ 0.74. Neglecting 
the variation of the bracket in the denominator of equation (|58D, we get for > 100 days 
F u (r ssa <C 1) oc £-°- 31 -°- 52e . The evolution of the optically thin flux is therefore considerably 
steeper if e = 1, -Fj,(r ssa <C 1) oc t~ 0,83 , compared to e = 0. Comparing with the observations 
after ~ 100 days there is much better agreement with the n re \ oc p W md V 2 model, compared 
to the n re i oc p model. The n re i oc p w i n d model gives too slow a decline for reasonable 



values of pi and m. This is consistent with the conclusions in section |4.3| . We also note the 
insensitivity of the optically thin flux to the magnetic field when p\ m 2, F v oc B^ 1 " 2 '' 2 ^ 
explaining the difficulty of separating the B oc i?" 1 and B oc t^ 1 cases. 

The break in the light curves at ~ 100 days occurs in our models because of the 
decrease in velocity, rather than as a result of a transition to the adiabatic phase. The 
rather abrupt change in the light curves, as well as in n ie \ in figure |TT|, shows that the 
transition from constant velocity to a more rapid decrease is fairly sharp, and should reflect 
a similarly well-defined flattening in the ejecta structure. 

For our parameters adiabatic cooling starts to dominate at t ~ 1000 days. However, as 
can be seen from figure O, the spectral index even at 3000 days is considerably higher than 
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the input spectral index, p [n = 2.1. In the adiabatic phase the optically thin flux scales as 

F u (T ssa « 1) OC r (Pi-l)/2-(l+2 £ )(l-m)_ (5g) 

We have also experimented with different slopes of the injection spectrum, p x . The 
most sensitive part of the spectrum is not surprisingly the epoch later than ~ 100 days, 
when most of the spectrum is optically thin. At these epochs the light curves are quite 
sensitive to pi, and variations in the range 2.0 ;$ p\ ^ 2.2 change the fit to the light curves 
markedly (as shown by the xl)- Our conclusion is therefore that the spectral index of the 
injected electrons should be in the above range, and not be changing in time. A value 
Pi = 2.1 is remarkably close to what diffusive shock acceleration gives, while it is not obvious 
why other acceleration mechanisms, like turbulent acceleration behind the shock, should 
give the same result. 

Another conclusion is that we can check the validity of the approximate equation ([43]) 
for the electron spectrum compared to the full solution of equation fl40|). We find that the 
difference between these approaches is small in both the magnetic field and non-thermal 
electron density. This is illustrated by the fact that our approximate analytical solutions 
(0) and Q51f ) accurately describe the numerical light curves. The results in section f| should 
therefore not be affected by the additional effects included in this section. 

6. FREE-FREE ABSORPTION, SYNCHROTRON SELF- ABSORPTION, 
AND THE RAZIN EFFECT FOR RADIO SUPERNOVAE 

One may ask under what circumstances free-free absorption, synchrotron self- 
absorption, or the Razin effect will be important. This is discussed extensively from the 
point of view of different supernova categories by Chevalier (1998). Here we limit ourselves 
to a few supplemental remarks. The radial optical depth to free-free absorption is 

(^f 2 ( i<1 ; hr , f (dy" < 6o > 

In the synchrotron cooling domain, and assuming for simplicity that p\ — 2, the synchrotron 
optical depth is r s « 1.6 x 10 4 / rel M (B 10 /30 G) 1 / 2 (^/2 x 10 4 km s" 1 )" 1 (t/10 days)~ 5/2 A 7/2 . 
From this we note that r s in this case is fairly insensitive to the magnetic 
field. In the adiabatic regime we get r s w 1.9 x 10 5 f ve \ M (Sio/30 G) 2 {V/2 x 
10 4 km s -1 ) -1 (t/10 days)~ 3 A 3 . In accordance with our results in section [4.2| , we 
further assume that the magnetic field is some fraction of the equipartition value, 



-27- 



B 10 » 1 x 10 2 (f) l l 2 M. G. The ratio of the synchrotron self-absorption and free-free optical 
depths is then in the cooling case 

^^6.7xlO 3 / r ci0 1/4 ^- 3/4 f^^y /2 ( £ r V f ^ V 2 A 3/2 , (61) 

r ff U0 5 K7 V2 x 10 4 km s" 1 / \ 10 days/ ' V ; 

and in the adiabatic case 



•l 



^ „ 4.8 x 10= /„, 4, (^f ( 2xi ; tms _, l A. (62) 



With the same assumptions, the ratio of the synchrotron cooling time to the expansion 
time is ^ 

^ « 1.1 x 1(T 2 0- 3/4 M- z ' 2 { -4 ) A 1 / 2 . (63) 



Unless M. <C 1 and/or ^<1, synchrotron cooling is likely to be important. 

Equation (|61| ) shows that in the cooling case 1) free-free absorption decreases in 
importance with time relative to synchrotron self-absorption, and 2) that, while insensitive 
to the magnetic field, the ratio is highly sensitive to the velocity and temperature of the 
circumstellar medium. 

The magnetic field is probably related to the progenitor system, and may differ between 
normal Type II, Type lib and lb supernovae. Because the latter types are likely to be 
components in binary systems, this fact may influence the circumstellar magnetic field. In 
addition, the turbulent amplification behind the shock may depend on the properties of the 
shock (e.g., whether it is cooling or not). Because of the lack of knowledge of this quantity 
it is fortunate that the optical depth is not very sensitive to the magnetic field. 

The circumstellar temperature is especially sensitive to the expansion velocity, since 
this, together with the mass loss rate, determines the shock temperature, and therefore also 
the EUV and soft X-ray flux from the supernova (Lundqvist & Fransson 1988). A higher 
expansion velocity leads to a higher shock temperature and a stronger UV and soft X-ray 
spectrum from Comptonized radiation, and thereby to a higher circumstellar temperature. 
Observationally, both SN 1993J and SN 19941 (Filippenko et al. 1995) had higher expansion 
velocity ( > 2 x 10 4 km s" 1 ) than SN 1979C and SN 1980K (~ 1.2 x 10 4 km s" 1 ). In 
addition, the color temperature of the initial burst of the supernova is higher for a compact 
progenitor, like for a Type Ib/c or lib supernova, which increases the heating of the 
circumstellar medium further. The difference in the supernova properties is illustrated 
in the different circumstellar temperatures found for SN 1979C and 1980K, compared to 
SN 1993J. In the former cases Lundqvist & Fransson (1988) found a peak temperature of 
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~ 2 x 10 5 K, while in the latter case the peak temperature was <^ 10 6 K (FLC96). Even 
more important, at 100 days the temperature was still ~ 2 x 10 5 K for SN 1993J, while it 
was only (2 - 3) x 10 4 K for SN 1979C and SN 1980K. 

Therefore, equations (fOD and (|62"D show that a higher expansion velocity increases the 
importance of synchrotron self-absorption relative to free-free absorption both directly and 
through the temperature. 

One should note that even if synchrotron self-absorption is larger than free-free, the 
rise of the light curve may be dominated by the free- free, since this acts externally to the 
synchrotron emitting region. Because ts itself is so sensitive to the velocity and temperature 
(equation jSOj) , the sensitivity to the velocity and progenitor structure is still present. The 
relative importance of synchrotron self-absorption and free-free absorption may therefore 
mainly be a result of the higher expansion velocity for the Type lib and Ib/c supernovae 
compared to ordinary Type II supernovae. 

For objects which are dominated by free-free absorption, equation floT| ) gives rise to 
the question whether synchrotron self-absorption may be important for the underlying 
spectrum also in these sources. Because of the limited number of observed frequencies this 
question is difficult to answer, although the above estimates makes this plausible. 

That a high expansion velocity is indicated for the supernovae showing synchrotron 
self-absorption has earlier been argued by Shklovskii (1985) and more recently by Chevalier 
(1998). These authors have discussed especially the radio emission from the Type lb 
supernovae, and conclude that the expansion velocities of these indeed have to be very high. 
The sensitivity to the temperature has, however, been discussed little before. 

As was illustrated in section |3.1.1| , the main problem for the Razin effect to be 
important is that it requires a very low magnetic field. To give a detectable radio emission 
it then follows that the column density of relativistic particles has to be unrealistically high. 
The high thermal plasma density needed for the Razin effect to dominate also means that 
Coulomb losses will become important. We therefore do not expect the Razin effect to be 
very important for supernovae showing an observable radio emission in other than quite 
exceptional cases. 



7. GENERAL CONSEQUENCES FOR VLBI OBSERVATIONS 



Both the dominance of the synchrotron self-absorption and the importance of electron 
energy losses have consequences for the VLBI observations. At wavelengths A ^> A Tssa= i the 
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intensity will be that of a uniform disk, while at A <C A Tssa= i limb brightening is important 
(Marscher 1985). For this latter situation the full set of equations (|I])-(0) should be used. 
Free-free absorption with the consequent limb darkening is only important for the longest 
wavelengths. VLBI observations by Bartel et al. (1994) and Marcaide et al. (1997) have 
been obtained at 1.35 and 2.01 cm from day 30 to day 91, at 3.56 cm from day 50 to day 
427, and at 6.16 cm from day 91 to day 1304. We find that the synchrotron optical depth is 
larger than unity up to day 16, 32, 63, and 125 at 1.35, 2.01, 3.56 and 6.16 cm, respectively. 
At times earlier than these epochs a uniform disk is therefore to be preferred as input to 
the VLBI modeling, while a limb brightened shell is more applicable after these epochs. 

Synchrotron losses have been shown to be important throughout the whole period of 
observations. If acceleration takes place close to the circumstellar shock in a very thin 
region, the emitting region will have a wavelength dependent thickness 

AR « oc VB-^X 1 ' 2 (64) 

4 

In this ideal case the thickness would therefore increase with wavelength by a factor of 
~ 40 % from 3 cm to 6 cm. One would also predict that the thickness would increase with 
time as the synchrotron timescale increases (figure 0). If the acceleration process takes place 
in the turbulent region behind the shock, possibly connected to the contact discontinuity, 
a more uniform emissivity would result. As long as synchrotron cooling dominates, the 
emissivity is not sensitive to the magnetic field distribution, because the flux is instead 
determined by the available energy in relativistic electrons. The opposite is true if the 
electrons are adiabatic. 

For the wavelengths where the supernova is optically thick to synchrotron self- 
absorption, one expects the polarization to be low, while in the optically thin phase the 
polarization can be high, unless the magnetic field is highly random. Because of the 
high degree of spherical symmetry the integrated polarization is likely to be low. VLBI 
polarization maps would therefore be highly interesting, although difficult to obtain in 
practice. A few resolution elements over the disk would suffice for this purpose. 

8. CONCLUSIONS 

Here we summarize our main conclusions: 

1. The light curves are characterized by synchrotron self-absorption in combination with 
external free-free absorption. While giving an excellent fit to the spectrum, the Razin 
effect can be excluded, based on an unrealistic injection efficiency 
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2. The circumstellar density is consistent with a constant mass loss rate and wind 
velocity, p wind oc r" 2 . 

3. The free- free absorption inferred from the combined synchrotron-self absorption plus 
free-free model agrees well with that estimated previously by FLC96. 

4. The magnetic field in the synchrotron emitting region can be reliably determined, and 
decreases like B oc R^ 1 , or alternatively as B oc The strength argues strongly 
for turbulent amplification behind the shock. The energy density is somewhat weaker 
than equipartition, u-q j 'i^herm ~ 0.14. 

5. Both synchrotron and Coulomb losses are important for the integrated electron 
spectrum, which cannot be characterized by a simple power law. 

6. The injected electron spectrum has a power law slope pi = 2.1, indicating that 
diffusive shock acceleration is operating. 

7. The number density of relativistic electrons scales as a fixed fraction of the thermal 
energy density behind the shock. This may argue for an injection mechanism which 
is controlled by the available thermal energy behind the shock. While the energy 
density in the relativistic electrons is only ~5x 10~ 4 of the thermal energy density, 
the ions may have a total energy density comparable to equipartition. 

8. The light curves, as well as the evolution of the injection rate of relativistic electrons, 
reflect the dynamics of the shock wave and therefore also the ejecta structure. 
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Fig. 1. — Angle averaged synchrotron spectra, including the Razin effect, for different power 
law indices, a. The Razin wavelength is given by Ar — 1-56 x 10 9 B/n c cm. 
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Fig. 2. — Observed VLA spectra at four epochs, from Van Dyk et al. (1994), together 
with fits based on a Razin suppressed synchrotron spectrum (solid line), external free- free 
absorption (dashed line), and synchrotron self-absorption (dotted line). 
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Fig. 3. — Derived values of Aro cc B/n c and A(t) from the fits to the Razin model. The 
solid line gives a best linear fit to the log _E>/n e versus logt relation. 
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Fig. 4. — Evolution of A Tssa= i and the scaled source function S'(t) = 7rS(t) R(t) 2 s /D 2 for the 
synchrotron self-absorption model with external free-free absorption. 
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Fig. 5. — Evolution of the magnetic field and column density of relativistic electrons in the 
synchrotron self-absorption plus external free-free absorption model. 
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Fig. 6. — Evolution of the mean intensity, J, at the shock, for -Eb-v — 0.08, M = 
5 x 1CT 5 M yr' 1 , and u w = 10 km s _1 . The quantity 16n 2 R^J is essentially the total 
luminosity at the shock. 
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Fig. 7. — Lorentz factor, 7 sy nch ) 7comp, 7coui, for which the synchrotron, inverse Compton, 
and Coulomb loss times, respectively, are equal to the expansion time. The dot-dashed 
line shows the region below which Coulomb losses are more important than synchrotron 
losses. Below the line A - B - C Coulomb losses dominate, and above the line A - B - D 
inverse synchrotron losses dominate, while the region C - B - D is adiabatic. The typical 7 
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Fig. 8. — Fits to the VLA spectra at six different epochs, including all suppression processes, 
as well as a self-consistent calculation of the electron spectrum. The dashed line shows 
the intrinsic spectrum from the shock, while the solid line also include external free-free 
absorption. 
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Fig. 9. — The free-free optical depth at 1 cm determined from the fits in this paper, compared 
to the optical depth from the model calculations of the circumstellar temperature in FLC96 
(dashed line). 
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Fig. 10. — Evolution of the magnetic field as a function of shock radius. The dashed line 
shows the expected evolution if B 2 oc PwmdV^ 2 , while the dotted line shows the evolution of 
an B oc Rr 1 field. 
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Fig. 11. — The number of relativistic particles as a function of the shock radius. The dashed 
line shows the expected evolution if n Te \ oc p W mdK, 2 , while the dotted line shows the evolution 
if n re i oc p wind . 




Fig. 12. — Calculated and observed radio light curves for the model with pi = 2.1, e = 1, 
/ re i = 1.85 x 10~ 4 and n = 6.5, discussed in the text, together with VLA observations, and 
the 2 cm Ryle observations (filled triangles). 
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Fig. 13. — Integrated line-of-sight electron spectra at 10 (solid), 100 (dotted), and 1000 
(dashed) days. The flattening of the spectrum caused by Coulomb losses below 7 ~ 10 is 
evident at especially 10 and 100 days, while above this energy synchrotron losses prevail. At 
10 days synchrotron self-absorption heating gives a weak bump in the spectrum at 7 ~ 10, 
corresponding to the energy of the electrons at which the optical depth is larger than one. 
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Fig. 14. — The effective spectral index p e g of the electron distribution at the Lorentz factor 
(7 = 85(-BA) -1 / 2 ) responsible for the emission at 1.3 - 21 cm for the model in figure 12. 



